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Abstract : Linear rate equations are used to describe the cascading decay of an initial 
| heavy cluster into fragments. This representation is based upon a triangular matrix of 



transition rates. We expand the state vector of mass multiplicities, which describes the 

process, into the biorthonormal basis of eigenmodes provided by the triangular matrix. 

When the transition rates have a scaling property in terms of mass ratios at binary frag- 

^ | mentation vertices, we obtain solvable models with explicit mathematical properties for 
O ' 

CNj ■ the eigenmodes. A suitable continuous limit provides an interpolation between the solvable 

t-H \ 

models. It gives a general relationship between the decay products and the elementary 

On 

transition rates. 

In a previous paper 1 -* we considered binary fragmentation processes where any frag- 
ment with mass number k breaks into fragments with mass numbers j and k — j, j = 
l,2.../c — 1, with a probability Wjk per unit of time. It was assumed that Wjk was time 
independent. By definition, wjk = if j > k and Wjk is symmetric if j is replaced by k — j, 
naturally. (For technical reasons, the coefficient Wj^j is twice the actual transition rate.) 
Let Nj (t) be the multiplicity of fragment j at time t in a process initiated from the decay 
of a cluster A, namely iVj(O) = 5jA- The model under study is described by the following 
set of linear, first order differential equations, 

alN A ' 7_1 w 

_i = _ CjNj + w jk N k , j = 1, ...A, Cj =J2^Y , (!) 

k=j+l £=1 

With components Nj, j = 1, ...A, for a column vector \Af >, the system, Eqs.(l), boils 
down to d\M > jolt = W\J\f > with a triangular matrix W. The solution of Eqs.(l) is 
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obviously \Af >= exp(— tW)\No >, a sum of exponentials whose rates of decay in time are 
the trivial eigenvalues of the triangular W, namely the diagonal matrix elements — c\. The 
purpose of the present letter is to take advantage of the expansion of the evolution on the 
biorthonormal basis of eigenmodes of W, 

A A 

\N >=J2 ' A > e M-tcx) < A|A/"o >= l A > exp(-c A t) Yj£, (2) 

A=l A=l 

where Y% is the last component of the bra eigenvector < A|. We are specially interested in 
the analytical structure of the bra eigenvectors < A|, because of the synthetic information 
carried by scalar products M\ =< X\Af >. 

In our previous investigations 1 - 1 we found numerical evidence for specific exponents 
g(A) which govern a power-like increase of the components of the bra eigenvectors, Y- oc 
ji( x ) . We took advantage of the fact that W has a fixed bra eigenvector < 1 |, whose 
components Yj = j, j = 1, ...A, express mass conservation. Indeed the total mass is 
Mi =< 1 \Af >= Y,f =1 jNj, with dMx/dt =< 1 \W\Af >= 0, since the first eigenvalue 
—c\ identically vanishes. It was also convenient to define the "mass weighted multiplicity" 
(MWM) vector U with components Ui = iNi, whose evolution is governed by a matrix 
V, with matrix elements Vjk = jV^jk/k, hence dU/dt = VU. A numerical analysis of U 
revealed scaling properties of the fragmentation process. The present paper attempts to 
give analytical proofs where only numerical evidence was previously found. In particular 
we are interested in those mathematical properties of the bra < A| eigenvectors which may 
enlighten the comprehension of the expansion, Eq.(2). 

Our main tools in this analysis are the following equations for the components Y { x of 
the bra eigenvectors of W, 

J-2 

Y t x w tj + Y^ wj-u = ( Cj - c x ) F/, j > A . (3) 

i=\ 

Similar equations govern the components of the ket eigenvectors, naturally. We keep 
in mind that Y^ = if j < A and X x = if i > A. A suitable (biortho) normalization 
results from the following boundary conditions, 

Y x x = X x = 1. (4) 
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Note that the components of the biorthonormal eigenvectors of the matrix V which governs 
the MWM are easily deduced, as = XY^/j and = iX^/X. 

It happens frequently that Cj is a monotonically increasing function of j. The physical 
interpretation is that the heavier the system, the larger its total decay rate. Under such a 
condition, the coefficients Cj —c\ in the right-hand side of Eqs.(3) are always positive when 
j > X. The left-hand sides of the same, Eqs.(3), depends only on positive coefficients Wij 
and components of rank i lower than j, X < i < j — 1. It is easy to calculate Y£ +1 from 
y A A alone, which was chosen to be positive from the normalization Y£ = 1. Thus Y^ +1 is 
positive. A trivial recursive argument then shows that all non vanishing components Y^ 
are positive. 1 

In the same way as in our earlier 1 -* studies we first select values of Wij which express 
scaling properties at vertices of binary fragmentation. For this, the family of models we 
consider 2 ^ uses the form Wij = [f(i/j)+f(l—i/j)] j~ a , and is parametrized by a "splitting" 
function / and an overall exponent a. For most of the present analytical study, we restrict 
a to the value a = 1, because it allows a convenient continuous limit, valid for finite values 
of A and large sizes A of the matrix. The value a = 1 is for instance the proper exponent 
for QCD fragmentation 3 ' 4 ^ when we consider energy-momentum fragmentation rather than 
mass fragmentation. 

A simplified, albeit generic form of / is parametrized by a second exponent, 6, in the 

form 



The variable x = i/j (with its mandatory symmetry complement 1 — x) is obviously the 
scaling mass ratio at the fragmentation vertex. We note that the cases a = 1, 6 = 
and a = 1, b = — 1 generate the same eigenvectors, since the corresponding matrices are 
strictly proportional, namely Wij\b=o = 2 Wij\b=-i = 

1 As a consequence of the positivity of the bras and of biorthogonality, the ket eigen- 
vectors must show highly oscillatory components. This demands a further, specific study. 
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For a = 1, b = — 1, the simple nature of Wij, namely Wij = reduces Eqs.(3) into 



1 ( Ci _! - c x ) Y x _ x + %1 = ( Cj - c x ) Y x . (3a) 



Hence a recursion relation for the bra eigenvector components is found, 

Y x i + X-1 
Y^~ i-X ' 

and accordingly the bra eigenvector components are given by the formula, 



(6) 



(i-A)! (2A-1)! K ) 

For large values of j, the recursion relation expands in powers of when j is large, 
according to 

% = 1 + ^ + , (8) 

% j 

which means that Y t x oc i 2X ~ l . This proves analytically that the scalar products < X\jV >= 
Yl,t=x Y^Ni of the bra eigenvectors with M express moments, of order q(X) = 2X — 1, of 
the multiplicity distributions, restricted to fragments at least as heavy as A. This justifies 
analytically, for a = 1 and b = —1,0, the numerical evidence obtained in our previous 
paper 1 ). 

For a = 1, b = 1 the rates become = j — The matrix V is simpler than W, 
since then Vy = — i) in the upper-right triangle, which does not vanish. The diagonal 
is unchanged, and so is, naturally, the lower- left, vanishing triangle. A straightforward but 
slightly tedious argument shows that 

One deduces from Eq.(9) that 



Y Ai _ U + 1) _ , , - 



y t - W+ l-A)- 1 + 7 + ° (1/i ) ' (10) 

hence now q(X) = X. This makes a second case where the power law, Y x oc j q ^ x \ can be 
proved analytically. 



If a is frozen to the value a = 1, it is tempting to interpolate between the behavior 
q(X) = 2A — 1 found for 6 = — 1 and that one, q(X) = A, valid for 6=1. Because of the 
identity of eigenvectors for the 6 = — 1 and 6 = parabolic interpolation 

A(17-46' 2 )+46 /2 -9 1 ( ,,, 

9int(A) = , 6=6+-, (11) 

is found. 

This interpolation turns out to be consistent with the continuous limit condition 
considered earlier, see Eq.(9) in reference 1 ), 

i=i i=i (12) 

= J / [x- 6 + (l-x)- b ] = / cte [1 - - (1 - z) 9 ] 

2 7o Jo 

where we take advantage of the symmetry Wij = Wj-ij. We can use this continuous 
limit, Eq.(12), to estimate q as a function of A when c\ is known, whether analytically or 
numerically. The identity of the cases 6 = and 6 = — 1 for Eq.(12) is transparent. Also, 
for any value of 6, the pair A = l,q = 1 solves obviously Eq.(12). Then, for 6 = 0, one 
gets at once from Eq.(12) the condition (A — 1)/A = 1 — 2/(^ + 1), hence q(X) = 2A — 1, 
in agreement with Eq.(8). Furthermore for 6=1, with the ansatz q(X) = A, the recursive 
relation, derived from Eq.(12), 

1 f 1 , x x + (1 - x) x - x x+1 - (1 - x) x+1 
c A+ i c\ = — = / dx 



X Jo x 
is easy to verify. 

We show on Fig.l the plots (full lines) of the values of (?(2), q(5) and q(8), respectively, 
obtained numerically from Eq.(12) for —1.5 < 6 < 1.5. The same Fig.l shows (dotted lines) 
the corresponding interpolation predictions qi n t(X) from Eq.(ll). As long as —1.2 < 6 < 
1.2, the continuous limit results and the parabolic interpolations do not differ by more 
than 4%. Significant deviations between such estimates occur, however, when 1.2 < |6|. It 
may be pointed out that the continuous limit estimates are expected to be more rigorous. 



Finally on Fig. 2 we show, for e.g. b = —0.5, b = 0.5 and b = 1.5, that, when provided 
by Eq.(12), the behavior of q as a function of A remains essentially linear. Namely a 
parametrization q(X) ~ a(b)X + (3(b) is reasonable, despite the fact that a(b) and (3(b) are 
more complicated than the quadratic forms present in Eq.(ll). 

In conclusion, there is a firm ground for the validity of the power law which governs 
the bra eigenvector components and the corresponding interpretation in terms of moments 
of the multiplicity vector Af. The apparent simplicity of the empirical law q(X) ~ a(b)X + 
/3(b), see for instance Eq.(ll), expresses a non trivial relationship between very different 
quantities, namely experimentally observable moments on one hand and, on the other hand, 
much less observable parameters b of the microscopic mechanism of binary fragmentation. 
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Figure Captions 

Fig.l : Behavior of the "eigenexponents" q(X) for the 2nd, 5th and 8th bra eigenvectors, 
as functions of the parameter b (model with scaling properties for fragmentation vertices) . 
Full lines: continuous limit predictions. Dotted lines: quadratic interpolations between 
the exact values, known analytically for b = —1, 0, 1. 

Fig. 2 : Strongly linear behavior of q(X) as a function of A, for different values of b. Full 
lines are drawn to guide the eye between integer values of A. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9412026vl 



This figure "figl-2.png" is available in "png" format from: 
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